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ABSTRACT 

The  finite  and  infinite  horizon  time  probability  of  ruin  are  important 
parameters  in  the  study  of  actuarial  risk  theory.  This  paper  introduces 
procedures  for  directly  estimating  these  key  parameters  from  a  random  sample 
of  observations  without  assumptions  as  to  the  parametric  form  of  the 
distribution  from  which  the  observations  arise.  The  estimators  introduced 
apply  to  most  of  the  classical  models  in  which  ruin  probabilities  are  used  and 
also  apply  to  a  much  broader  class  of  models.  The  procedures  are  based  on  the 
concept  of  sample  reuse,  an  old  idea  in  statistics  which  is  becoming  more 
popular  due  to  the  widespread  availability  of  high  speed  computers.  In  this 
paper,  the  almost  sure  consistency  of  the  estimators  is  established.  Further, 
finite  sample  properties  of  the  estimators  are  investigated  in  a  simulation 
study . 


AMS  (MOS)  Subject  Classifications:  Primary  62P05;  Secondary  60G40,  62G99. 
Key  Words:  reverse  martingale,  bootstrap  , 


Work  Unit  #4  (Statistics  and  Probability) 


SIGNIFICANCE  AND  EXPLANATION 


The  probability  of  ruin  is  the  probability  that  claims  against  a  risk¬ 
taking  enterprise  exceed  its  initial  capital  or  reserve  plus  income  at  some 
point  of  time.  If  ruin  occurs  before  a  fixed  time,  such  as  one  or  ten  years, 
this  is  called  the  finite  horizon  time  problem  and  otherwise,  the  infinite 
horizon  time  problem.  Even  in  the  simplest  models,  ruin  probabilities  can  be 
difficult  to  compute.  There  is  a  vast  literature  on  this  aspect  of  the 
problem  dating  back  at  least  to  the  early  1900's. 

In  this  paper  it  is  shown  how  to  approximate  the  probability  of  ruin 
directly  from  available  claims  data.  Various  properties  of  the  approximation 
procedure  are  provided.  Despite  the  long  history  of  research  in  this  area, 
it  is  not  surprizing  that  these  procedures  have  not  been  developed  previously 
because  they  are  computer  intensive  and  would  have  been  computationally 
burdensome  ten  years  ago.  However,  in  today's  era  of  fast  computation,  these 
procedures  are  inexpensive  to  implement.  An  example  of  the  procedure  is 
provided  via  simulating  claims  data  and  calculating  approximations  to  the 
probability  of  ruin  of  the  enterprise. 
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NONPARAMETRIC  ESTIMATION  OF  THE  PROBABILITY  OF  RUIN 
Edward  W.  Frees 

§1.  INTRODUCTION 

Let  (X^,Y^ )  i  -  1,2,...  be  l.i.d.  random  vectors  with  joint  bivariate  distribution 

function  F.  Use  the  random  variable  to  represent  the  it*1  claim  amount  and  the 

nonnegative  random  variable  Y ^  to  represent  the  ith  interarrival  time  between  claims. 

For  t  >  0,  define  the  number  of  claims  by  time  t  as 

N( t )  -  l  I( Y,  +  .. .  +  Yk  <  t) 
k>  1 

where  I ( • )  is  the  indicator  function.  Premiums  are  assumed  to  arrive  at  a  known  steady 
rate,  say,  P  per  unit  time.  Thus,  the  amount  that  claims  exceed  premiums  by  time  t  is 

N(t) 

U( t )  -  l  xk  -  Pt  .  (1.1) 

k-1 

N(t) 

Here  we  interpret  the  sum  £  to  be  zero  when  N(t)  »  0.  We  are  primarily  interested  in 

k-1 

the  probability  that  U(t)  exceeds  an  initial  reserve  u  at  some  time  t  prior  to  or 

at  T,  the  horizon  time.  This  probability  can  also  be  defined  by 

<i>(U,T)  -  p(  sup  U(t)  >  u)  .  (1.2) 

0<t<T 

We  are  also  interested  in  the  probability  that  U(t)  eventually  exceeds  an  initial 
reserve  u, 

ip(u)  -  ( u ,<*» )  =  lim  i|>(u,T)  .  (1.3) 

>p->oo 

These  probabilities,  t|i(u,T)  and  if(u),  are  called,  respectively,  the  finite  and  infinite 
horizon  time  probability  of  ruin. 

The  probability  of  ruin  is  a  key  parameter  in  the  collective  theory  of  risk  and  has 
received  considerable  attention  over  the  years.  For  several  different  types  of 
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introduction  to  this  theory,  see  Feller  (1971),  Beard  et.  al.  (1984),  Gerber  (1979)  and 
Biihlmann  (1970).  The  behavior  of  the  finite  horizon  time  probability  of  ruin  has  been 
widely  investigated  not  only  because  of  its  practical  relevance  but  also  because  explicit 
cor.putat  ion  of  iHu,T)  is  difficult  except  in  the  most  trivial  cases.  Explicit 
calculations  have  been  given  for  some  specific  forms  of  F,  cf . ,  Sell  (1978)  and  Thorin 
and  wikstad  (1976).  However,  use  of  these  explicit  solutions  has  been  limited  due  to  their 
complexity  and  their  dependence  on  a  specific  form  of  F.  Because  of  these  difficulties, 
papers  giving  approximations  of  \)i(u,T)  suggested  by  limit  theorems  (e.g.,  as 
u  *  ■»,  T  ♦  oo)  are  abundant  in  the  literature.  Hie  most  successful  of  these 
approximations  seem  to  be  the  diffusion  approximations  given  by  Siegmund  (1979)  and  applied 
by  Asmussen  (1984).  See  Lalley  (1984)  for  a  refinement  of  Siegmund's  work.  Another  class 
of  methods  for  calculating  i)i(u,T)  is  the  Monte-Carlo  method.  Surprisingly,  the  Monte- 
Carlo  method  has  received  only  limited  attention  in  the  risk  theory  literature.  For  some 
accounts,  see  Beard  et.  al .  (1984)  and  Seal  (1978).  Note  that  this  method  does  not  depend 
on  a  specific  form  of  F  but  does  depend  on  complete  knowledge  of  F. 

The  approximations  of  ip (u ,T)  and  i)»(u)  given  in  this  paper  are  different  in  nature 
from  those  sketched  above  and  are  inspired  by  the  concept  of  sample  reuse.  Sample  reuse  is 
an  old  idea  in  statistics,  popularized  in  the  nineteen-forties  by  Hoeffding  (1948)  and  more 
recently  by  Efron  (cf.,  1982).  The  idea  for  our  applications  is  as  follows.  Consider  the 
random  variable,  Z  =  suPo<t<T  U<t)-  Since  the  distribution  of  U(t)  is  completely 
determined  by  F,  then  the  distribution  of  Z  is  completely  determined  by  F.  Note  that 
from  (1.2),  \Mu,T)  =  1  -  P(Z  <  u).  By  the  usual  multivariate  Glivenko-Cantelli  Theorem, 
knowledge  of  (Xi,Yi)  i  =  1,2,...  completely  determines  F  and  thus,  for  a  sufficiently 
large  sample  size  n,  from  { (X^ ,Y^  ) } "_ ^  we  can  build  a  good  approximation  to  F,  the 
usual  (multivariate)  empirical  distribution  function.  This  raises  the  natural  question  of 
how  to  build  a  reasonable  estimator  of  i|>(u,T)  based  on  a  random  sample  of  size  n, 

( (X  , Y  ) }"  .  The  purposes  of  this  paper  are  to  argue  that  this  is  an  important  question 
and  to  develop  estimators  of  i|/(u,T)  (and  \|/(u)  )  that  possess  desirable  properties. 


In  §2,  estimators  of  iMu,T)  and  iji(u)  are  introduced.  The  almost  sure  (a.s.) 
insistency  of  these  estimators  is  proved  in  §3.  A  small  simulation  study  is  given  in 
§4.  Because  of  the  nature  of  the  approximations,  the  assumptions  made  in  our  development 
are  different  and  more  general  than  the  usual  risk  theory  assumptions.  In  §5  we  discuss 
these  differences  and  make  other  concluding  remarks. 


§2.  ESTIMATORS 


Let  ( (* i t > ^  a  random  sample  from  a  population  with  distribution  function 

F.  Let  (aij  >a2j  *  •  •  •  'anj  )  =  A.j  be  the  jth  permutation  of  {l,2,...,n},  ]  =  1,...,nl  We 

intend  to  reuse  the  sample  by  considering  reordered  pairs  {(X^  (Y^  > } ^  r  j  *  1,...,n! 

aij  tj 

Define  sn  =  Y1  *  •••  +  yn‘  For  t^ie  3*"^  permutation,  the  number  of  claims  by  time 


NA,(t)  =  V  I(Y-  +  ...  +  Y,  <  t)  t  <  S„ 

3  kt,  1 j  aXj  " 


t  >  Sn  , 


and  the  amount  that  claims  exceeds  premiums  by  time  t  is 

NA  (t) 

3 

UA  j  ( t )  -  l  X  -  Pt  . 

J  k-1  k3 

For  the  jth  permutation,  we  have  ruin  if 

I (sup  UAj(tl  >  u)  ■  1  (2.1) 

where  the  supremum  is  over  the  set  {t  :  0  <  t  <  min(Sn,T)}.  Note  that  to  compute  the 

function  in  (2.1)  one  does  not  need  to  evaluate  UA^(t)  at  all  t  e  ( 0 ,min ( Sn ,T) )  but 

only  at  the  random  time  points  0,  Y  ,Y  +  Ya  , . . .  ,Y  +  ...  +  Y  .  The  firBt  type 

1j  a1j  2 j  °1j  nj 

of  estimator  we  consider  is  the  average  over  all  permutations, 

t|in<u,T)  =  (n!)-1  I  J  (sup  UAjlt)  >  u)  (2.2) 

and 

*n(u)  «  *n(u,Sn)  (2.3) 

where  \  is  the  sum  over  all  permutations  of  {1,2,...,n}.  Note  that  we  could 
alternatively  define 

t)in(u)  »  ♦  (u,*)  =  lim  i)i  (u,T)  . 

The  consistency  of  these  estimators  is  provided  in  the  following 


Theorem  2.1. 


Suppose  that  EY  >  0.  Then,  for  each  T, 

lim  i^lu.T)  »  i|i(u,T)  a.s.  (2.4) 

n** 

and 

lim  ^(u)  =  4i(u)  a.s.  (2.5) 

n-*» 

By  the  SLLN,  the  requirement  EY  >  0  ensures  that  Sn  ♦  00  a.s.  The  proof,  given  in 

§3,  is  based  on  the  idea  that  can  be  shown  to  be  a  reverse  martingale  plus  negligible 

terms.  We  remark  that  the  estimators  defined  in  (2.2)  and  (2.3)  each  require  the 

evaluation  oi  nl  indicators  of  ruin,  an  extensive  amount  of  computations  even  for 
moderate  sample  sizes  (say,  n  >  10).  Because  of  this  computational  difficulty,  we  now 
introduce  a  bootstrap  estimator  of  the  probability  of  ruin.  The  bootstrap  methodology,  a 
name  coined  by  and  a  methodology  popularized  by  Efron  (cf.,  1982),  is  also  computer¬ 
intensive  but  does  not  require  a  prohibitive  amount  of  computation. 

Let  B  «  B(n)  be  a  positive  integer  depending  on  n  such  that  B  +  “  as  n  + 

Based  on  the  observed  sample  ^  #  we  draw  B  bootstrap  realizations  of  <(ln  in 

the  following  two  steps.  For  b  «  1,...,E, 

Step  1.  Make  n  independent  (conditional  on  {(X  ,Y  )}”^)  draws  with  replacement 

from  {(X^.Y^)}"^  to  get  (X^.Y^)  for  i  «■  1,...,n. 

Step  2.  Define 

N*b(t)  “  l  I(Y*b  +  ...  +  Y*b  <  t) 
k“1  ' 

ar.d  compute 

*b 

•b  N  (t)  *b 

V  (u,T)  =  l(  sup  (  l  K  -  Pt)  >  u)  . 

0<t<min(S  ,T)  k= 1  ’ 

n 

The  bootstrap  estimates  are  defined  by 

-1  B  *b 

t*(u,T)  =  B  l  !(/  ( u ,  T)  (2.6) 

n  ,  .  n 

b»  1 


-5- 


a 


♦•(u)  *  »J>*(u,Sn)  .  (2.7) 

The  consistency  of  these  estimators  is  provided  in  the  following 
Theorem  2.2. 

Suppose  F'Y  >  0  and 

log  n  =  0(B1/2)  •  (2.8) 

Then,  tor  each  T, 

lim  4**(u,T)  =  ^(u,T)  a.s.  (2.9) 

n  x*5 

and 

lim  ^*(u)  =  4Mu)  a.s.  (2.10) 

n-  " 

The  proof  of  Theorem  2.2  is  given  in  §3.  The  condition  on  B  in  (2.8)  guarantees  that 
B  grows  sufficiently  quickly  to  achieve  a.s.  convergence.  In  §4  we  investigate  the  finite 


sample  properties  of  i)/*  is  a  small  simulation  study 


§3.  PROOFS 

Define  Y1n  >'f2n  '  *  •  •  '^nn  to  the  or<^er  statistics  of  Y.j  ,Y2  ,  •  •  •  ,Yn  and  let  X^n 
the  claims  amount  associated  with  Yln,  i  =  1,...,n.  Let  Gn  =  o((X^n,Y^n), 
i  -  (Xi(Y|),  i  >  n)  for  n  >  1,  a  non increas ing  sequence  of  sub  o-fields. 

Because  all  the  arguments  follow  easily  for  the  case  T  =  *,  we  only  give  proofs  for  the 
a.s.  consistency  of  =  ^n(u,T).  We  now  define  a  version  of  and  show  that  it  is  a 

reverse  martingale.  Later  we  show  that  this  version  is  close  to  i|<n  in  the  appropriate 
sense  . 

For  the  ]th  permutation  of  {1,2, . ..,n),  let 


NA j ( t )  =  NAjft)  +  l  I(Sk  <  t) 
k>n 

be  a  version  of  the  number  of  claims  by  time  t,  j  =  1,...,nl  With  a^j  *  k  for 
k  >  n,  define 

NA' (t) 

4>*  ( u.T)  =  (n!)_1  )  l(  sup  l  X  -  Pt  >  u)  .  (3.1) 

n  p  0<t<r  k»?  kj 

We  have  the  following  property  for  this  version  of  the  probability  of  ruin. 

Lemma  3 . 1 

For  each  u,T,  (t'(u,T),G  )  is  a  reverse  martingale, 
n  n 

Vroof  : 

It  is  easy  to  see  that  0»*  =  (u,T)  is  G  -measurable  and  G  -integrable .  With 

n  n  **  n 

C(r)  defined  as  in  (1.1),  we  have  the  relation 

=  e(i(  sup  U(t)  >  u)|G  }  .  (3.2) 

"  0<t<T 

Thus , 

eU^G  .)  =  E1EU(  9UP  U(t>  >  u)lGn^Gn+lt 

0<t<T 

*  e(i(  sup  U(t)  >  u)|G  ,(  -  >f;+1  .  + 

0<t<T 


be 
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i 

t 

Proof  of  Theorem  2.1: 

From  (3.2),  we  have  that  EyMu.T)  =  v(u,T).  From  Lemma  3.1  and  the  (reverse) 

martingale  converyence  theorem.  It  is  easy  to  show  that 

lira  i|i  1  ( u  ,T)  =  i(i(u,T)  a.s.  (3.3) 

n—  " 

Define  the  stopping  time 

n 

T  =  inf (n  >  1  -.  l  (X.  -  PY  )  >  u|  . 

k=1  *  * 

Now,  it  is  easy  to  see  that  4'n(u,T)  and  i(j^(u,T)  differ  only  on  the  set  {n  <  x  <  “}  . 

Thus 

0  <  i|j  (u,T)  -  1(1*  (u,T)  =  I(n  <  T  <  “)  ♦  0  a.s. 
n  n 

This  and  (3.3)  are  sufficient  for  the  proof.  If 
Proof  of  Theorem  2.2: 

Denote  =  iJi*(u,T).  From  (2.4),  a  sufficient  condition  for  (2.9)  is,  for  E  >  0, 

ln  P<  It'S  -  *nl  >£)<». 

We  show 

l  P<+*  -  >  e)  <  «  ,  (3.4) 

the  proof  of  the  other  inequality  being  similar. 

Define  Fn  =  offx^.y^),  i  =  1,...,n)  for  n  >  1,  a  nondecreasing  sequence  of  sub 
o-fields.  Now,  from  the  condition  (2.8),  there  exists  a  positive  constant  K  such  that, 
for  sufficiently  large  n, 

n'1  >  exp(-KB1/2)  .  (3.5) 

From  the  Markov  inequality,  with  6  >  G  \  we  have 

PCk*  -  \||  >  e)  <  exp{-e6KB1,/2}E{exp(6KB1/2(i)'*  -  ^  )  ) }  • 

n  n  n  n 

1  /2 

From  (3.5),  exp{-E|5KB  }  is  summable.  Thus,  to  prove  (3.4),  we  need  only  show 

sup  E(exp(  sB1^2  ( 1(1*  -  <*/))}<”  (3.6) 


■8- 


where  s  =  6k.  From  (2.6),  since  conditionally  on  Fn,  41*  is  the  mean  of  a  binomial 

random  variable,  we  have 

E{exp(sB1//2(t*  -  ))} 

n  n 

=  E((1  -  )exp{-s')i  B  ^2}  +  i(j  exp{s(1  -  41  )  B  1,2})B 
n  n  n  n 

-  E(1  +  s2*  (1  -  <|/  )/<2B)  +  0(b'3/2))B 
n  n 

<  (1  +  s2/(8B)  +  0(b'3/2))B 

by  a  Taylor-series  expansion.  This  is  sufficient  for  (3.6)  and  hence  (2.9).  ^ 


§4.  SIMULATION 


In  this  section,  finite  sample  properties  of  the  bootstrap  est imators  introduced  in  §2 
were  investigated.  A  simple  example  was  used  so  that  calculation  of  exact  probabilities  of 
ruin  and  compatison  with  other  studies  were  possible.  Claim  amounts  were  assumed  to  be 
exponential!/  distributed  with  mean  1.  The  claims  were  assumed  to  arrive  as  a  Poisson 
process  with  intensity  P  =  .8,  i.e.,  interarrival  times  are  independent  and  exponentially 

distributed  with  mean  1.25.  The  claims  amount  and  arrival  times  were  assumed  to  be 
independent  and  premiums  arrive  with  unit  intensity  (P  =  1).  Under  these  assumptions,  it 
can  be  shown  that 

Hi(u)  =  .8  exp{-.2u}  ,  (4.1) 

thus  giving  an  easy  expression  for  exact  values  of  p(u).  furthermore,  in  a  recent  study, 
Asmussen  (  1984)  has  provided  exact  values  of  p(u,T)  and  several  popular  approximations  of 
i|i(u,T)  for  various  values  of  T. 

The  bias  (BIAS)  and  root  mean  square  error  (RMSE)  were  used  to  judge  the  performance 
of  the  estimators.  All  computations  were  done  on  a  VAX  11/750  owned  and  operated  by  the 
Department  of  Statistics  at  the  University  of  Wisconsin-Madison.  The  IMSL  Fortran 
subroutines  produced  the  random  deviates. 

In  Table  1  we  give  the  results  of  the  performance  of  the  bootstrap  estimator  of  the 
infinite  horizon  time  probability  of  ruin.  The  tables  give  the  criteria  for  ruin 
probabilities  ip  ( u)  =  1%,  5%,  10%,  40%  and  for  sample  sizes  n  =  10,  30,  50,  100,  150. 

The  ruin  probabilities  were  chosen  to  represent  a  range  which  is  typically  of  interest  to 
the  actuarial  community.  The  sample  sizes  were  selected  to  represent  small  and  moderate 
numbers  of  claims.  For  small  sample  sizes  n  =  10,  30,  50,  B  =  200  bootstrap 
replications  were  run  to  compute  the  estimators  for  each  simulation  trial.  For  larger 
sample  sizes  n  =  100,  150,  we  used  B  =  100  bootstrap  replications  to  compute  the 
estimators.  Using  a  smaller  B  produced  considerable  savings  in  computer  run  time  and  did 
not  seem  to  affect  significantly  the  accuracy  of  the  results. 

In  Table  2  we  give  the  results  of  the  performance  of  the  bootstrap  estimator  of  the 
finite  horizon  time  probability  of  ruin.  The  number  of  bootstrap  replications  are  as  in 
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table  _ 


1>(u) 


1* 


5% 


10% 


40% 


n 


B 


BIAS 


10 

30 

50 

100 

150 


200 


100 


10 

30 

50 

100 

150 


200 


100 


-.00902 

-.00773 

-.00364 

.00470 

.01365 

-.04716 

-.02602 

-.01025 

.01065 

.02400 


10 

30 

50 

100 

150 


200 


100 


-.09023 

-.03965 

-.01650 

.00960 

.02000 


10 

30 

50 

100 

150 


200 


100 


-.24040 

-.06825 

-.04225 

-.04010 

-.02250 


RMSE 

,00999 

.01770 

.02116 

.03910 

.05615 

.04061 

.06063 

.07711 

.10527 

.11267 

.09648 

.10944 

.12108 

.14964 

.14489 

.32200 

.26922 

.24297 

.22975 

.21270 
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TABLE  2  -  bOOTbTRAP  ESTIMATOR  OF  »(u,T) 


T 

\F(u,T) 

n 

B 

BIAS 

RMSE 

13.8 

.00007 

30 

200 

.00118 

.00546 

50 

.00138 

.00459 

100 

100 

.00043 

.00346 

150 

.00048 

.00272 

200 

.00023 

.00173 

41.3 

.00145 

30 

200 

.00005 

.00581 

50 

.00330 

.01590 

100 

100 

.00805 

.02699 

150 

.00670 

.02276 

200 

.00370 

.01144 

68.8 

.00338 

30 

200 

-.00188 

.00611 

50 

.00140 

.01566 

100 

100 

.00892 

.03525 

150 

.01387 

.04527 

200 

.01017 

.03118 

96.4 

.00491 

30 

200 

-.00341 

.00674 

50 

.00013 

.01560 

100 

100 

.00739 

.03489 

150 

.01509 

.05172 

200 

.01579 

.04764 

As  was  expected,  the  performance  of  the  infinite  time  estimator  was  better  the  closer 
♦lu)  was  to  50%.  In  Table  1,  the  asymptotic  theory  comes  quickly  into  play  when 
iKu)  =  40%  as  evidenced  by  the  decreasing  BIAS  and  RMSE  terms  with  increasing  sample 
size.  For  i|/(u)  »  1o%,  we  see  some  levelling  off  of  the  KMSE  term  from  n  =  100  to 
n  =  150.  The  smaller  probabilities  showed  no  evidence  of  levelling  off  for  the  sample 
sizes  considered.  Although  the  motivation  of  decreasing  B  at  higher  sample  sizes  was  to 
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decrease  computing  costs,  this  did  not  seem  to  affect  the  performance  of  the  estimators. 

See  Efron  (1982)  for  further  discussion  of  the  selection  of  B. 

Perhaps  the  most  interesting  fact  of  the  simulation  study  was  that  the  performance  of 
the  finite  time  estimators  improves  as  T  decreases.  In  Table  2,  the  asymptotic  theory 
comes  quickly  into  play  when  T  m  13.8.  For  T  ■  41.3,  we  see  the  decreasing  trend  in  the 
BIAS  and  RMSE  terms  beginning  at  a  larger  sample  size  n  *  100.  The  intuition  is  as 
follows.  For  T  »  13.8,  on  the  average  it  requires  about  11  observations  to  check  for 
ruin  by  time  13.8  (since  11  *  1.25  *  1  “  13.75).  With  a  sample  size  of  n  *  100,  we  have 
approximately  9  independent  and  identical  realizations  of  an  indicator  of  ruin  by  time 
T  »  13.8.  Repeating  this  reasoning  for  T  “  41.3,  the  reader  can  verify  that  we  have  only 
approximately  3  independent  and  identical  realizations  of  this  indicator  of  ruin,  while 
the  goal  of  the  estimators  introduced  in  §2  is  to  reuse  sampling  information,  by  increasing 
the  horizon  time  T  we  increase  the  number  of  observations  necessary  to  check  for  ruin  and 
thus  reduce  the  number  available  for  resampling. 

55.  DISCUSSIOM 

The  estimators  introduced  in  52  are  defined  in  terms  of  the  classical  method  of 
collecting  data,  i.e.,  based  on  a  random  sample  of  size  n.  However,  it  is  easy  to  modify 
these  estimators  for  other  methods  of  data  collection  and  retain  their  important 
statistical  properties.  For  example,  suppose  an  insurance  company  would  like  to  estimate 
the  probability  of  ruin  based  on  one  year  of  observed  data.  Then  the  sample  size  itself  is 
random.  However,  it  is  easy  to  see  that  the  usual  results  on  random  change  of  time  (cf., 
Csorgo  and  Rlvesz,  1980,  Theorem  7.1.1)  can  be  applied  to  preserve  the  properties  stated  in 
Theorems  2.1  and  2.2.  In  the  risk  and  queuing  literature,  this  change  of  time  is  usually 
refered  to  as  a  transition  to  operational  time . 

The  most  important  drawback  of  the  estimators  i|»n  and  'I'*  is  that  they  rely  heavily 
on  the  independence  of  the  bivariate  pairs  (Xi,Yi).  While  this  assumption  is  used  in  most 
models  constructed  to  calculate  the  probability  of  ruin,  other  models  such  as  a  model  which 


■ 


uses  the  mixed  Poisson  process  for  the  claims  number  process  (of.,  beal  and  Gerber,  1984) 

1o  not  . 

Except  for  the  imports  c  assumption  of  independence,  the  estimators  and  1(1*  and 

'heir  jropert.ies  are  relatively  free  from  assumptions  when  compared  to  other  estimators  of 
the  probability  of  ruin.  Both  estimators  are  nonparametr ic  in  the  sense  that  they  do  not 
assume  a  particular  parametric  form  nor  knowledge  of  the  distribution  function  F.  In  the 
assumptions  we  have  not  precluded  the  case  P(Y^  =  0)  >  0,  thus  allowing  for  the 
possibility  of  multiple  claims  at  any  point  in  time.  Further,  we  have  not  assumed 
independence  between  the  claim  arrival  time  and  the  claim  amount.  If  this  assumption  is 
nade ,  we  conjecture  that  a  different  estimation  procedure  can  be  constructed  that  uses  the 
iata  more  efficiently,  in  some  sense.  We  leave  this  as  an  open  question  for  future 
research . 

We  remark  here  that  questions  of  finding  the  optimal  or  most  efficient  procedure,  in 
some  sense,  have  not  been  addressed.  For  example,  in  Step  1  of  the  definition  of  4»*  in 
§2,  if  the  drawing  of  { was  made  without  replacement  then  inspection  of  the 
proof  of  Theorem  2.2  shows  that  Theorem  2.2  still  holds.  Whether  the  drawing  should  be 
made  with  or  without  replacement  we  leave  as  an  area  of  future  research.  The  bootstrap 
estimator  i()*  is  computationally  similar  to  the  usual  Monte-Carlo  procedure.  However,  it 
is  different  philosophically  in  that  the  Monte-Carlo  procedure  assumes  knowledge  of  the 
underlying  distribution  function  while  the  bootstrap  does  not.  This  similarity  suggests 
that  dynamic  factors  in  ruin  probabilities  such  as  interest,  inflation,  economic  cycles, 
etc.,  that  have  been  incorporated  in  Monte-Carlo  methods  ( cf . ,  Beard  et.  al.,  1984,  Chapter 
7)  may  be  incorporated  in  bootstrap  estimators.  We  leave  this  as  an  open  area  for  future 
research . 
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